Generalidades

En esta entrega se describe la estrategia de calibración usada con 5 estaciones de caudal. Los datos observados de caudal han sido previamente tratados (ver: Tratamiento de caudales).

Periodo de calibración y validación

La siguiente tabla muestra los periodos de validación.

t23 = c("2004-01-01","2005-12-31") 
t30 = c("2005-01-01","2006-12-31") 
t36 = c("2008-01-01","2009-12-31") 
t38 = c("2007-01-01","2008-12-31") 
t41 = c("2006-01-01","2007-12-31")
T.valida = rbind(t23,t30,t36,t38,t41)
kable(T.valida)

Las siguientes figuras muestran en gris los periodos de validación, el resto de la serie se utiliza para calibración.

#lectura de datos
qobs <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/qobs_5est_2004_2012.RDS")

# ajuste de caudal
qobs$q23 = 1.630*qobs$q23
qobs$q30 = 0.789*qobs$q30
qobs$q36 = 1.184*qobs$q36
qobs$q38 = 0.724*qobs$q38
qobs$q41 = 1.120*qobs$q41

#conversion a mm
qobs$q23 = convertFlow(qobs$q23, from="cumecs", to="mm", area.km2 = 687)
qobs$q30 = convertFlow(qobs$q30, from="cumecs", to="mm", area.km2 = 1092)
qobs$q36 = convertFlow(qobs$q36, from="cumecs", to="mm", area.km2 = 2744)
qobs$q38 = convertFlow(qobs$q38, from="cumecs", to="mm", area.km2 = 3159)
qobs$q41 = convertFlow(qobs$q41, from="cumecs", to="mm", area.km2 = 4896)

fig = plotly_empty()

fig = fig %>% add_segments(x = ~t23[1], y=max(qobs$q23,na.rm = T), yend = max(qobs$q23,na.rm = T), 
                           xend = t23[2], name="Validación Paso Pache",
                           line = list(color = 'black', width = 2))

fig = fig %>% add_segments(x = ~t30[1], y=max(qobs$q30,na.rm = T), yend = max(qobs$q30,na.rm = T), 
                           xend = t30[2], name="Validación Paso Roldán",
                           line = list(color = 'blue', width = 2))

fig = fig %>% add_segments(x = ~t36[1], y=max(qobs$q36,na.rm = T), yend = max(qobs$q36,na.rm = T), 
                           xend = t36[2], name="Validación Fray Marcos",
                           line = list(color = 'red', width = 2))

fig = fig %>% add_segments(x = ~t38[1], y=max(qobs$q38,na.rm = T), yend = max(qobs$q38,na.rm = T),
                           xend = t38[2], name="Validación San Ramon",
                           line = list(color = 'orange', width = 2))

fig = fig %>% add_segments(x = ~t41[1], y=max(qobs$q41,na.rm = T), yend = max(qobs$q41,na.rm = T),
                           xend = t41[2], name="Validación Paso Pache",
                           line = list(color = 'green', width = 2))

fig <- fig %>%add_lines(x = ~qobs$Rdate, y = ~qobs$q23, type="scatter",mode = "lines",
               name="Paso de los Troncos", line = list(color = 'black', width = 2))

  
fig <- fig %>%add_lines(x = ~qobs$Rdate, y = ~qobs$q30, type="scatter",mode = "lines",
               name="Paso Roldán", line = list(color = 'blue', width = 2))

  
fig <- fig %>%add_lines(x = ~qobs$Rdate, y = ~qobs$q36, type="scatter",mode = "lines",
               name="Fray Marcos", line = list(color = 'red', width = 2))

  
fig <- fig %>%add_lines(x = ~qobs$Rdate, y = ~qobs$q38, type="scatter",mode = "lines",
               name="Paso San Ramón", line = list(color = 'orange', width = 2))

  
fig <- fig %>%add_lines(x = ~qobs$Rdate, y = ~qobs$q41, type="scatter",mode = "lines",
               name="Paso Pache", line = list(color = 'green', width = 2))



fig <- fig %>% layout(#legend = list(x = 1, y = 0.9),
                      yaxis = list(title="Runoff (mm)",
                                   showticklabels = TRUE, ticks="outside",
                                   showline=TRUE, showgrid=TRUE),
                      xaxis = list(title="",showticklabels = TRUE, 
                                   ticks="outside",showline=TRUE, 
                                   showgrid=TRUE, range=as.Date(c("2004-01-01", "2012-12-31"))),
                      height=300,
                      width=700)

  
fig